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Abstract. We present oxygen abundances for 72 F and G dwarf stars in the solar neighbourhood. Using the kinematics of 
the stars we divide them into two sub-samples with space velocities that are typical for the thick and thin disks, respectively. 
The metallicities of the stars range from [Fe/H] » -0.9 to +0.4 and we use the derived oxygen abundances of the stars to: (1) 
perform a differential study of the oxygen trends in the thin and the thick disk; (2) to follow the trend of oxygen in the thin disk 
to the highest metallicities. We analyze the forbidden oxygen lines at 6300 A and 6363 A as well as the (NLTE afflicted) triplet 
lines around 7774 A. For the forbidden line at 6300 A we have spectra of very high S IN (>400) and resolution (R > 215 000). 
This has enabled a very accurate modeling of the oxygen line and the blending Ni lines. The high internal accuracy in our 
determination of the oxygen abundances from this line is reflected in the very tight trends we find for oxygen relative to iron. 
From these abundances we are able to draw the following major conclusions: (i) That the [O/Fe] trend at super-solar [Fe/H] 
continues downward which is in concordance with models of Galactic chemical evolution. This is not seen in previous studies 
as it has not been possible to take the blending Ni lines in the forbidden oxygen line at 6300 A properly into account; (ii) That 
the oxygen trends in the thin and the thick disks are distinctly different. This confirms and extends previous studies of the 
other a-elements; (iii) That oxygen does not follow Mg at super-solar metallicities; (iv) We also provide an empirical NLTE 
correction for the infrared O i triplet that could be used for dwarf star spectra with a S IN such that only the triplet lines can be 
analyzed well, e.g. stars at large distances; (v) Finally, we find that Gratton et al. (1999) overestimate the NLTE corrections for 
the permitted oxygen triplet lines at ~ 7774 A for the parameter space that our stars span. 

Key words. Stars: fundamental parameters - Stars: abundances - Galaxy: disk - Galaxy: formation - Galaxy: abundances - 
Galaxy: kinematics and dynamics 



1. Introduction 

Oxygen is, next to hydrogen and helium, the most abundant el- 
ement in the Universe. Of the three stable isotopes lfi O, n O, 
and ls O, the first is the dominating one, making up ~ 99.8 % of 
the total oxygen content in the Solar system. Oxygen is a bona 
fide primary element that essentially only forms in the interiors 
of massive stars through hydrostatic burning of mainly He, C, 
and Ne. By analyzing elemental abundances in the atmospheres 
of long-lived F and G dwarf stars it is not only possible to de- 
termine the chemical composition of the gas that the stars were 
born out of but also to trace the chemical history in the different 
stellar populations in the Milky Way. In this respect the oxygen 
fossil record is of extra importance and is often used in models 
of Galactic evolution. It can, among other things, be used to 
measure the rates of supernovae type II (SN II) and supernovae 
type la (SNIa) with time (e.g. Wheeler at al. 1989). An over- 
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abundance of oxygen ([O/Fe] >0) indicates that the region 
have had a high star formation rate and undergone a fast chemi- 
cal enrichment (e.g. Tinsley 1979 , Matteucci & Greggio 1986 1. 

The determination of oxygen abundances is unfortunately 
often troublesome due to the limited number of available oxy- 
gen lines in the visual part of a stellar spectrum. The main 
indicators are the Oi triplet at ~ 7774 A, the forbidden [Oi] 2 
lines, especially those at 6300 A and 6363 A, the ultraviolet 
(UV), and the infrared (IR) OH lines. The analysis of these 
lines all have their difficulties. The triplet at 7774 A should 
be ideal to work with since its lines are strong and are lo- 
cated in a clean part of the spectrum that is free from blending 
lines. The lines are, however, strongly affected by deviations 
from local thermal equilibrium (LTE) (see e.g. Kiselman 1991 1 
and, even if these deviations are considered not very well un- 
derstood even for the Sun, they are probably due to convec- 
tive inhomogeneities (Kiselman 1993). The forbidden lines 
are very robust indicators but are hard to work with since 



1 Abundances expressed within brackets are as usual relative to so- 
lar values: [O/Fe] = log(N /N F X - log(N /N Fc ) Q 

2 These brackets indicate that the spectral line is forbidden and 
should not be confused with the notation of abundances given rela- 
tive to solar values. 
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they are both weak (the measured equivalent widths in the 
Sun are W^(6300)«5mA and W^(6363) * 3 mA, Moore et 
al. [T966). Both lines are also blended by lines from other 
species. The [Oi] line at 6300 A has two Nil lines in its right 
wing (Lambert ll978l Johansson et al. 2003 1 and the [Oi] line 
at 6363 A has a possible CN contribution (Lambert 1978 1. 

The conclusion from the Joint Discussion 8 during the IAU 
General Assembly in Manchester in 2000 (New Astronomy 
Reviews 45, 2001), which was devoted to a discussion of oxy- 
gen abundances in old stars, is that the most reliable indica- 
tors of stellar oxygen abundances are the forbidden [O i] lines 
(Barbuy et al. 200fl. The [Oi] line at 6300 A is blended and 
is located in a part of the spectrum that is severely affected by 
telluric lines. If we want to use this line for abundance analy- 
sis it is important to be able to take these effects into account 
properly. This is best done by using spectra with high signal- 
to-noise ratios (S /N) and high resolutions (R) together with ac- 
curate atomic line data, log gf- values in particular. 

In the Milky Way it is mainly the halo stellar population 
that shows the large overabundances of oxygen relative to iron, 
indicative of a fast star formation. There is however an ongoing 
controversy about the size of this overabundance. Depending 
on which indicators that are used in the determination of the 
oxygen abundances one can find different [O/Fe] trends for 
these stars. The forbidden oxygen line at 6300 A and the IR 
OH lines generally produce a constant value of [O/Fe] as 0.5 
for halo stars with [Fe/H] <-l, see e.g. Nissen et al. (|7002i, 
Gratton & Ortolani dl986> . Sneden et al. dl979> . and Lambert 
et al. ( 1974) for abundances from the [Oi] line at 6300 A, and 
e.g. Melendez & Barbuy (J2002 1, Melendez et al. (2001 1, and 
Balachandran et al. (120011 for abundances from the IR OH 
lines. Other indicators such as the Oi triplet at 7774 A and the 
UV OH lines give higher [O/Fe] values with an uprising trend 
when going to lower [Fe/H], see e.g. Israelian et al. (2001 ) and 
Abia & Rebolo ( 1989) for abundances from the Oi triplet, and 
Boesgaard et al. 0999) and Israelian et al. {1998 1 for abun- 
dances from the UV OH lines. However, agreement between 
abundances from the UV OH lines and the forbidden lines has 
also been found (e.g. Bessel et al. 119911 ) and also between the 
Oi triplet and the forbidden lines (e.g Carretta et al. 2000 
Mishenina et al. l2000b . 

For stars more metal-rich than the halo (i.e. [Fe/H]>-1 
to -0.5) the observed trends are further complicated by the 
fact that there is a mixture of stars from different stellar 
populations, each with their own chemical history (see e.g. 
Gratton et al. 2003). The observed [O/Fe] trend for the halo 
stars generally extends to [Fe/H] — 0.5 with an overabun- 
dance of [O/Fe] -0.4-0.5. Nissen & Schuster STWH . how- 
ever, found that the most metal-rich halo stars to be less en- 
hanced in [O/Fe] than the thick disk stars at the same metallic- 
ities. Quite good agreement between the different oxygen in- 
dicators are generally found at these metallicities (see works 
cited above). For the thick disk no (previous) study has in- 
cluded thick disk stars with [Fe/H] > -0.35 when studying 
the [O/Fe] vs. [Fe/H] trend. Below this metallicity, thick disk 
stars are generally found to show an [O/Fe] trend similar to 
the halo trend, i.e. a constant value of [O/Fe] ~ 0.4-0.5 (see 
e.g. Nissen et al. 120021 Tautvaisiene et al. 120011 Prochaska 



et al. l2UU01 Gratton et al. 2000 Chen et al. 2000). The thin 
disk stars present small to moderate overabundances of oxygen 
(0 < [O/Fe] < 0.2) for metallicities down to [Fe/H] « -0.8 with 
the largest values of [O/Fe] at the lower [Fe/H] values (see e.g. 
Reddy et al. 120351 Nissen et al. 12001 Edvardsson et al. 179951 
Nissen & Edvardsson 1992; Kjsergaard et al. 1982). The dis- 
tinction between the thin and thick disk is however somewhat 
unclear. While e.g. Prochaska et al. (2000) found them to be 
distinct in terms of abundances, Chen et al. ( 2000 ) found the 
elemental abundance trends in the thin and thick disks to fol- 
low smoothly upon each other. 

At super-solar metallicities ([Fe/H] > 0) a constant value 
of [O/Fe] ~ was found by Nissen & Edvardsson ( 1992), and 
Nissen et al. (2002). Feltzing & Gustafsson ( 1998 ) has a some- 
what larger scatter in the [O/Fe] trend and also find a potential 
decline towards higher [Fe/H]. A decline in [O/Fe] is supported 
by the studies of Castro et al. J1997I and Chen et al. ( 2003 ), but 
is at variance with other studies discussed above. 

In the solar neighbourhood it is also possible to come across 
stars which might originate in the Galactic bulge or inner disk 
(Pompeia et al. 2002). They are generally old (> lOGyr) and 
Pompeia et al. ( 2002 ) found that they show an overabundance 
of oxygen when compared to their disk counterparts for metal- 
licities ranging -0.8 < [Fe/H] < +0.4. 

In this paper we will determine the oxygen abundances for 
thin and thick disk stars with -0.9 < [Fe/H] < +0.4. The stars 
have been selected on purely kinematical grounds (see Bensby 
et al. 2003 for a full discussion of the selection process) in or- 
der to address two important questions: 1) The [O/Fe] trends 
at [Fe/H] < in the thin and thick disks; 2) The [O/Fe] trend at 
super-solar metallicities. 

Our oxygen abundances are based on the forbidden [Oi] 
line at 6300 A. In order to properly account for the Ni 
blends and the fact that the line generally is very weak we 
have obtained spectra of very high quality (/?>215 000 and 
S /N >4QQ). We have used the new log gf- values for the Ni 
blends from Johansson et al. ( 2003 ) who also showed that the 
Ni blend actually consists of two isotopic Ni components. 

The paper is organized as follows. In Sect. [2] we describe 
the observations and the data reductions. In Sect.|5Jwe describe 
the atomic data used in the abundance determination, which in 
turn is described in Sect. Section [5] deals with the errors in 
the resulting abundances. The abundance trends that we derive 
from the different indicators are presented in Sect. [6] In Sect.Q 
we put our results into the contexts of Galactic chemical evo- 
lution. We finally give a summary in Sect. [8] 

2. Observations and data reductions 

2.1. Stellar sample 

The stellar samples have been kinematically selected to contain 
stars with high probabilities of belonging to either the thin or 
the thick Galactic disks. These probabilities were derived un- 
der the assumption that the stellar space velocities Ulsr, Vlsr, 
and Wlsr in the thin and thick disks have Gaussian distribu- 
tions. For each star it is then possible to calculate the probabil- 
ity that it kinematically belongs to either the thin or the thick 
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Fig. 1. a) The CES image slicer profile, b) A zoom in on the 
straylight pedestal beneath the light pattern of the image slicer 
profile. 



disk. Since we are working with nearby stars we also need to 
take the local number density of the thin and thick disk stars 
into consideration. For the thick disk we used a density of 6 % 
in the solar neighbourhood and consequently 94 % for the thin 
disk. We then formed the TD/D ratio (i.e. the thick disk prob- 
ability divided by the thin disk probability) for each star. All 
thin disk stars in the sample (apart from three) have thin disk 
probabilities that are more than ten times larger than their thick 
disk probability (i.e. TD/D < 0.1). Three thin disk stars (Hip 
3170, 95447, 100412) have higher TD/D ratios (but still ap- 
proximately two times more probable of belonging to the thin 
disk than the thick disk). The majority of the thick disk stars 
are more than ten times more probable to belong to the thick 
disk than to the thin disk (i.e. TD/D > 10). For a few thick disk 
stars we loosened this criterion somewhat in order to increase 
the statistics, allowing thick disk stars with 1 <TD/D< 10 in 
the sample. Since the two thick disk sub-samples (TD/D> 10 
and 1 < TD/D < 10) show exactly the same trends for other el- 
ements such as Na, Mg, Al, Si, Ca, Ti, Cr, Ni, and Zn (see 
Bensby et al. 2003), we will not distinguish between them in 
this study. 

Apart from the three thin disk stars with slightly larger 
TD/D ratios, two additional thin disk stars (Hip 5862 and Hip 
15131), and one additional thick disk star (Hip 15510) the stel- 
lar sample is the same as in Bensby et al. (2003) (i.e., 66 stars 
in common). The sample now consists of 22 thick disk stars 
and 50 thin disk stars (see Table [Q. The selection criteria and 
the stellar samples are fully described in Bensby et al. (J2003). 
The additional 6 stars are described in a forthcoming paper 
(Bensby et al., in prep.) where we moreover analyze a further 
30 thin and thick disk dwarf stars in the northern hemisphere, 
observed with the SOFIN spectrograph on the NOT telescope 
on La Palma. 
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Fig. 2. Example of the interference pattern that originates from 
reflections within the optical fiber. It is best seen in a featureless 
spectrum, in this case a fast rotating B star. The two telluric 
lines are the same as those in Fig. [3] The spectrum has been 
wavelength calibrated but not normalized and not corrected for 
the radial velocity of the star. 

2.2. Observations 

Observations were carried out with ESO's 3.6 m telescope 
on La Silla in Chile during two observing runs in September 

2000 (4 nights, TB and SF as observers) and August/September 

2001 (6 nights, TB as observer). We used the Coude Echelle 
Spectrograph (CES) at its highest resolution, R a 215 000, cen- 
tered on the forbidden oxygen line at 6300 A (hereafter denoted 
as [Oi]63oo)- This setting gives a spectrum with a wavelength 
coverage of ~40A. All spectra have a S /N of at least 400. 
Long exposures times (given in Tabled were usually split into 
three exposures. Since telluric lines are present in this wave- 
length region we also observed rapidly rotating B stars. These 
spectra were then used to divide out telluric lines from the ob- 
ject spectra. 

Spectra covering the whole optical wavelength region from 
3560 A to 9200 A were also obtained with the FEROS spectro- 
graph on the ESO 1.5 m telescope during the same observing 
runs. These spectra have R ~ 48 000 and S/N> 150 and pro- 
vide Fe and Ni abundances as well as the atmospheric parame- 
ters for the stars (analysis presented in Bensby et al. 2003). 

Spectra for the three additional stars (Hip 5862, 15131, 
15510) were obtained with the UVES spectrograph on the 
VLT/UT2 telescope Kueyen on Paranal in Chile in July 

2002 (TB as observer). These spectra have R~ 110000 and 
S/N ~ 350 (Bensby et al., in prep.). 

2.3. Reduction of CES spectra 

The reduction procedure from raw images to wavelength cali- 
brated spectra was carried out with the MIDAS 3 software. The 
reductions are straightforward using standard routines concern- 
ing cosmic ray filtering, background subtraction, flatfield di- 
vision, wavelength calibration, and continuum fitting. All ex- 

ESO-MIDAS is the acronym for the European Southern 
Observatory Munich Image Data Analysis System which is developed 
and maintained by the European Southern Observatory. 



4 



T. Bensby et al.: Oxygen trends in the Galactic thin and thick disks 



0.8 



0.8 



0.6 



a) 


1 1 


i i i | i i i i 






[oi] \ 1 - 


c 




-H 2 


1 1 1 1 II 1 1 1 II 1 1 1 

b) 






[01] \ / - 






i i 1 i , J i 



I I I i I I I I i i I I I Lj i I 

6298 6299 6300 6301 



Wavelength [A] 

Fig. 3. Example of the removal of lines originating in Earth's 
atmosphere. The telluric lines, O2 at 6298.5 A, and O2+H2O at 
6299.2 A are indicated as well as the [O i]63oo line, a) Spectrum 
of Hip 75181 prior to removal of telluric lines and b) after re- 
moval of telluric lines by division with a spectrum of a fast 
rotating B star. The spectrum has been wavelength calibrated 
but not corrected for the radial velocity of the star. Note the 
telluric oxygen emission line at ~ 6300.25 A. 



posures for a single star were reduced individually and then 
the extracted spectra were co-added before the setting of the 
continuum. For those stars that had a telluric absorption line 
coinciding with the [Oi]63oo line we divided the spectra with 
a featureless spectrum from a fast rotating B star. Finally the 
co-added spectra were shifted to rest wavelengths by measur- 
ing the doppler shift of eight evenly distributed Fe 1 lines that 
have accurate wavelengths from Nave et al. (119941 . There are, 
however, some properties of the CES spectra that are worth dis- 
cussing in further detail. 

Straylight pedestal: The background level is different on the 
two sides of the spectrum in the direction perpendicular to 
the dispersion direction (see Fig. [3- This is most likely dif- 
fuse light produced by the image sheer itself (Kiirster 2001 1. 
In the wavelength direction this background is slightly non- 
uniform. Actually it is a low-resolution version of the stellar 
spectrum that for all practical purposes can be treated as fea- 
tureless (Kiirster 2001 1. The [Oi]63oo line profile from individ- 
ual slices were checked to see if this straylight pedestal had any 
influence on the line profile and strength from different slices. 
We especially checked for variations in the line strength (equiv- 
alent width) of the [O1J6300 line. No discernible effects (apart 
from differences in the S /N) were, however, found between the 
slices. 



Interference pattern: A regular sinusoidal pattern can also 
be seen in the continuum of a featureless CES spectrum (see 
Fig. |2}. Most likely this has its origin from reflections within 
the optical fiber that feeds the spectrograph (Kiirster 2001 1. 
For low metallicity stars with a low density of spectral lines 
it is possible to perform a Fourier analysis of the pattern and 
remove it. However, most of our stars are rich in spectral lines 
which makes a removal nearly impossible. Our B star spectra 
indicate that the amplitude of this pattern is very small, typi- 
cally 0.2% of the continuum flux (see Fig. |2j. Unless spectra 
with S /N> 500 are striven for, it is therefore safe to ignore this 
feature in the reduction process. 

Telluric absorption lines: Telluric lines were divided out using 
spectra of fast rotating B stars with the IRAF 4 task telluric 
which is part of the noao package. The procedure takes care 
of small possible shifts in the dispersion zero-points as well 
as of the variations in the abundance of the telluric species by 
shifting and scaling the B star spectrum. An example of the 
removal of the telluric absorption lines is given in Fig. [3] Of 
the 72 stars in our sample there were eight stars (Hip 3086, 
3185, 3704, 15510, 75181, 82588, 88622, and 95447) that had 
a telluric absorption line that blended the [Oi]63oo line. Only 
for these eight stars did we divide out the telluric lines. 

Telluric emission lines: Unfortunately there were also occa- 
sions when the telluric oxygen emission line (see Fig. [3} co- 
incides with the stellar [Oi]63oo line. Since the presence and 
strength of this feature vary strongly with time it is very diffi- 
cult to remove. A possible way to remove it would be to ob- 
serve a part of the sky next to the object immediately after the 
observation of the star, and then use this sky spectrum to di- 
vide out the sky emission. We did, however, not obtain such 
sky spectra. For seven stars (Hip 3170, 9085, 12611, 22263, 
86731, 90485, 83229) we were unable to determine the oxy- 
gen abundance based on the [Oi]63oo line due to this telluric 
emission feature. 

3. Atomic data 

The oscillator strength of the weak Ni line at 6300.35 A has 
recently been measured in the laboratory by Johansson et 
al. (|2003 1. They showed that it actually contains two isotopic 
components, /l( 58 Ni) = 6300.335 A and ,l( 60 Ni) = 6300.355 A. 
Both components have log gf = — 2.11. Assuming a solar iso- 
topic ratio of 0.38 for the abundances of 6() Ni and 58 Ni and that 
these two isotopes make up 94.4% of the total Ni abundance 
(CRC 2003 1, the weighted log g/- values for these isotopes are 
logg/( 58 Ni) = -2.275 and log g/( 60 Ni) = -2.695. 

The log gf- values for the oxygen lines have been taken 
from the VALD database (Piskunov et al. 119951 Kupka et 
al. [19991 Ryabchikova et al. 119991 . The source for the [Oi] 

4 IRAF is distributed by National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation, USA. 
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Table 1. The 72 program stars. T e ff, \ogg, £ t , [Fe/H], and [Ni/H] are taken from Bensby et al. (2003). All abundances have been 
normalized to the Sun. Both NLTE corrected, according to the prescription in Gratton et al. ( 2000 1, and un-corrected abundances 
are given for the triplet lines. For oxygen abundances from the individual triplet lines see Table @] The last three columns give 
the spectrographs and exposure times used: C=CES (ESQ 3.6 m), F=FEROS (ESQ 1 .5 m), and U=UVES (ESQ VLT/UT2). 
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0.17 


0.23 


C(60) 


F(6) 


80337 


147513 


6094 


5880 


4.49 


1.10 


2.09 


0.03 


-0.04 


0.03 




-0.04 


C(20) 


F(3) 


80686 


147584 


6098 


6090 


4.45 


1.01 


2.28 


-0.06 


-0.13 


-0.05 


-0.09 


0.01 


C(20) 


F(2) 


81520 


149612 




5680 


4.53 


0.65 


1.70 


-0.48 


-0.49 






-0.24 




F(ll) 


83601 


154417 


6349 


6167 


4.48 


1.21 


3.27 


0.09 


0.00 


0.07 




0.04 


C(34) 


F(5) 


84551 


156098 


6409 


6475 


3.79 


2.00 


3.98 


0.12 


0.08 


0.12 




0.28 


C(18) 


F(3.3) 


84636 


156365 




5820 


3.91 


1.30 


2.68 


0.23 


0.32 


0.15 




0.28 


C(55) 


F(9) 


85007 


157466 




6050 


4.37 


1.10 


2.43 


-0.39 


-0.45 


-0.26 




-0.27 


C(75) 


F(10) 


85042 


157347 


6465 


5720 


4.40 


1.00 


1.65 


0.03 


0.06 


0.00 


0.06 


-0.05 


C(40) 


F(6.7) 


86731 


161239 


6608 


5840 


3.79 


1.43 


3.37 


0.25 


0.33 






0.26 


C(39) 


F(5) 


86796 


160691 


6585 


5800 


4.30 


1.05 


2.20 


0.32 


0.41 


0.14 




0.24 


C(27) 


F (2.5) 


87523 


162396 


6649 


6070 


4.07 


1.36 


2.17 


-0.40 


-0.39 


-0.31 


-0.31 


-0.23 


C(35) 


F(6) 


90485 


169830 


6907 


6339 


4.05 


1.55 


2.59 


0.12 


0.11 




0.06 


0.25 


C(31) 


F(5) 


91438 


172051 


6998 


5580 


4.42 


0.55 


1.62 


-0.24 


-0.26 


-0.18 




-0.23 


C(45) 


F(4) 


94645 


1 79949 


7291 


6200 


4.35 


1.20 


3.90 


0.16 


0.1 1 


0.06 


0.1 1 


0.23 


C (40) 




95447 


182572 


7373 


5600 


4.13 


1.10 


2.14 


0.37 


0.47 


0.20 


0.23 


0.38 


C(18) 


F (2.5) 


96536 


184985 


7454 


6397 


4.06 


1 .65 


2.72 


0.03 


0.00 


0.05 




0.20 


C(21) 


F(3.3) 


98785 


190009 


7658 


6430 


3.97 


1.90 


4.25 


0.03 


-0.03 


0.11 


0.01 


0.18 


C(60) 


F(7.5) 


99240 


190248 


7665 


5585 


4.26 


0.98 


1.72 


0.37 


0.44 


0.19 


0.17 


0.32 


C(4) 


F(l) 


100412 


193307 


7766 


5960 


4.06 


1.20 


2.18 


-0.32 


-0.36 


-0.18 




-0.03 


C(45) 


F(6) 


102264 


197214 




5570 


4.37 


0.60 


1.81 


-0.22 


-0.26 


-0.12 




-0.15 


C(80) 


F(ll) 


103682 


199960 


8041 


5940 


4.26 


1.25 


2.34 


0.27 


0.35 


0.14 


0.12 


0.21 


C (70) 


F (6.5) 


105858 


203608 


8181 


6067 


4.27 


1.17 


2.04 


-0.73 


-0.74 


-0.46 




-0.34 


C(30) 


F(l) 


109378 


210277 




5500 


4.30 


0.78 


1.71 


0.22 


0.26 


0.16 


0.16 


0.24 


C(60) 


F(8) 


110341 


211976 


8514 


6500 


4.29 


1.70 


3.63 


-0.17 


-0.24 


-0.11 


-0.12 


-0.01 


C(45) 


F(6) 


113137 


216437 


8701 


5800 


4.10 


1.16 


2.17 


0.22 


0.27 


0.13 


0.12 


0.27 


C(35) 


F (5.5) 


113357 


217014 


8729 


5789 


4.34 


1.00 


1.91 


0.20 


0.26 


0.10 


0.24 


0.11 


C(21) 


F(3.2) 


113421 


217107 


8734 


5620 


4.29 


0.97 


1.75 


0.35 


0.42 


0.16 




0.27 


C(45) 


F(6) 


117880 


224022 


9046 


6100 


4.21 


1.30 


3.16 


0.12 


0.14 


0.05 




0.16 


C (60) 


F(5.5) 



lines is the compilation of Wiese et al. (1966 1 who based 
the transition probabilities on theoretical calculations from 
Stoffregen & Derblom (fl960l, O mholt dl960l . Naqvi (IT95TV 
and Yamanouchi & Horie |1952J. For the Oi triplet lines the 
original sources are Biemont & Zeippen ( 1992) and Hibbert et 
al. (1T99H . 

The broadening of atomic lines by radiation damping was 
considered in the determination of the abundances and the 
damping constants (y ra d) for the different lines were collected 
from the VALD database. 

Collisional broadening, or Van der Waals broadening, was 
also considered. The width cross-sections are taken from 
Anstee & O'Mara lfl995l . Barklem & O'Mara (fT997l PT998 1. 



and Barklem et al. ( 1998 2000). Lines present in these works 
have been marked by an "S" in Table 12 For spectral lines not 
present in these works (marked by a "U" in Table |3 we apply 
the correction term (5yc) to the classical Unsold approximation 
of the Van der Waals damping. Following Mackle et al. d!975> 
this parameter was set to 2.5 for the lines considered here. 
The atomic data are summarized in Table [2] 

4. Abundance determination 

Oxygen abundances have been determined by line synthesis 
of the [Oi]63oo line (Sect. l4~4l i. by equivalent width measure- 
ments of the [Oi]6363 line (Sect. l4~5T > and the Oi triplet lines at 
~ 7774 A (Sect. 14. 6> . The actual generation of synthetic spectra 
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Identifications 






log« 


ft 


"tOX- macro 


[Fe/H] 


[Ni/H] 


[O/H] 


[O/H] 


[O/H] 


Spectrograph (EXPT1ME) 


Hip 


HD 


HR 


IK-1 


[cgs] 


[km s~'] 


[kms-'j 






6300 


6363 


7774 




[min] 














THICK DISK STARS 












3086 


3628 




5840 


4.15 


1.15 


1.90 


-0.11 


-0.12 


0.09 




-0.04 


C (100) 


F (10) 


3185 


3795 


1 73 


5320 


3.78 


0.70 


1.64 


-0.59 


-0.57 


-0.14 


-0.21 


-0.11 


C (45) 


F (6.7) 


3497 


4308 




5636 


4.30 


0.80 


1 .76 


—0.33 


—0.31 


—0.05 


-0. 10 


0.03 


C (45) 


F (10) 


3704 


4597 




6040 


4.30 


1 .08 


9 T() 


-0.38 


-0.40 


-0.08 




-0.25 


C (180) 


F (10) 


5315 


6734 




5030 


3.46 


0.86 


1.66 


-0.42 


-0.36 


0.00 


-0.01 


0.03 


C (45) 


F (10) 


14086 


18907 


914 


5 1 10 


3.5 1 


0.87 


1 .64 


—0.59 


-0.53 


—0.16 




-0.20 


C (42) 


F (5) 


15510 


20794 


1008 


5480 


4.43 


0.75 


1.28 


-0.41 


-0.36 


-0.01 


-0.10 




C (50) 


U (0.2) 


17147 


22879 




5920 


4.33 


1.20 


1 .74 


-0.84 


-0.83 


-0.32 




-0.32 


C (107) 


F(10) 


75181 


136352 


5699 


5650 


4.30 


0.78 


1.69 


-0.34 


-0.34 


0.00 




0.01 


C (22) 


F(5) 


79137 


145148 


6014 


4900 


3.62 


0.60 


1.53 


0.30 


0.39 


0.20 


0.13 


0.36 


C (28) 


F (4) 


82588 


152391 




5470 


4.55 


0.90 


2.53 


-0.02 


-0.06 


0.01 


0.00 


-0.01 


C (54) 


F(8) 


83229 


153075 




5770 


4.17 


0.97 


1.52 


-0.57 


-0.51 




-0.23 


-0.10 




F(ll) 


84905 


157089 




5830 


4.06 


1 .20 


1.98 


-0.57 


-0.56 


—0.15 




—0 17 


C(75) 


F(ll) 


88622 


165401 




5720 


4.35 


0.80 


2.10 


-0.46 


-0.45 


-0.05 


-0.09 


0.02 


C (60) 


F(9) 


96124 


183877 




5590 


4.37 


0.78 


1.78 


-0.20 


-0.17 


0.08 


0.08 


0.07 


C(80) 


F(13) 


98767 


190360 


7670 


5490 


4.23 


0.66 


1.98 


0.25 


0.31 


0.21 


0.13 


0.38 


C(25) 


F(5) 


103458 


199288 




5780 


4.30 


0.90 


1.69 


-0.65 


-0.64 


-0.29 




-0.29 


C (60) 


F(8) 


108736 


208988 




5890 


4.24 


1.05 


1.97 


-0.38 


-0.35 


-0.01 


-0.01 


0.06 


C(80) 


F(14) 


109450 


210483 




5830 


4.18 


1.10 


1.43 


-0.13 


-0.15 


0.05 


0.08 


0.06 




F(10) 


109821 


210918 


8477 


5800 


4.29 


1.05 


1.72 


-0.08 


-0.10 


0.10 




0.01 


C(45) 


F (6.7) 


110512 


212231 




5770 


4.15 


1.05 


2.12 


-0.30 


-0.28 


0.00 




0.07 


C(180) 


F(10) 


118115 


224383 




5800 


4.30 


1.00 


1.20 


-0.01 


-0.04 


0.04 




-0.01 




F(10) 



Table 2. Atomic line data for the different oxygen lines and 
the blending Ni lines. Column 1 gives the element, col. 2 
the wavelength, col. 3 the lower excitation potential, and 
col. 4 the correction factor to the classical Unsold damping 
constant. Column 5 indicates if the broadening by collisions 
have been taken from Anstee & O'Mara ( 1995 1, Barklem & 
O'Mara fl997l [19981 . and Barklem et al. ( fl998l l2000l l (indi- 
cated by an "S") instead of the classical Unsold broadening 
(indicated by an "U"). Column 6 gives the radiation damping 
constant and col. 7 the log gf- values (references are given in 
Sect.0. 



El. 


A 

(A) 


X\ 

(eV) 




DMP 


Trad 

(s" 1 ) 


logg/ 


[Ol] 


6300.304 


0.00 


2.50 


U 


1.0e+05 


-9.819 


[Oi] 


6363.776 


0.02 


2.50 


u 


1.0e+05 


-10.303 


Oi 


7771.944 


9.14 


2.50 


s 


4.8e+07 


0.37 


Oi 


7774.166 


9.14 


2.50 


s 


4.7e+07 


0.22 


Oi 


7775.388 


9.14 


2.50 


s 


4.5e+07 


0.00 


58 Ni 


6300.335 


4.27 


2.50 


u 


2.7e+08 


-2.275 


60 Ni 


6300.355 


4.27 


2.50 


u 


2.7e+08 


-2.695 



was done with the Uppsala Spectrum synthesis program and for 
the calculation of elemental abundances from measured equiv- 
alent widths we used the Uppsala Eqwidth abundance pro- 
gram. The stellar model atmospheres are standard 1-D LTE and 
were calculated with the Uppsala MARCS code (Gustafsson et 
al. [19751 Edvardsson et al. fT^l Asplund et al. nWfl . 

4.1. Stellar atmosphere parameters 

Stellar model atmosphere parameters for 66 of the 72 stars in 
the sample were derived in Bensby et al. (2003 1 using ~ 150 
Fe i lines measured in the FEROS spectra. The parameters for 
the additional 6 stars that are included here have been deter- 



mined in the same way (using UVES spectra for 3 of them) 
and will be presented in Bensby et al. (in prep.) (together with 
the SOFIN stars, see Sect. |3. In summary, the tuning of the 
model atmospheres includes the following steps: effective tem- 
peratures were determined by requiring Fe i lines with different 
lower excitation potentials to produce equal abundances; sur- 
face gravities were determined using Hipparcos parallaxes and 
stellar masses; microturbulences were determined by forcing 
the Fe i lines to give the same abundances regardless of line 
strength. The atmospheric parameters ([Fe/H], T e g, logg, and 
£ t ) for the whole sample are given together with the derived 
oxygen abundances in Tabled 

4.2. Macroscopic line broadening 

Apart from line broadening due to atomic properties the line 
profile is affected by the instrument profile, the line-of-sight 
component of the stellar rotation (v ro t ■ sin i), and the macrotur- 
bulence. The instrument profile broadening is set by the resolu- 
tion of the spectra (R). To determine the contributions from the 
other two broadening mechanisms we used the two strong Fe i 
lines at 6301.501 A and 6302.494 A (the former can be seen to 
the right in Fig.[3J. So, in order to broaden the spectral lines we 
first convolved the synthetic spectra with a Gaussian profile in 
order to produce the instrumental broadening and then with an 
elliptical profile to produce the combined effect of the stellar 
rotation and macroturbulence. Since we have chosen to include 
the macroturbulence broadening into the elliptical broadening 
of the rotation we refer to this profile as Vrot-macro- 

In Fig. 0] we illustrate the process for Hip 75181. The first 
step was to tune the log ^/-values of the two Fe i lines in order 
to make the shape of the synthetic line profiles similar to the ob- 
served lines. This was done since the Fe abundances that were 
used in the models have been determined from ~ 150 other, 
generally weaker, Fe i lines (see Bensby et al. 2003 1. We cre- 
ated 10 synthetic spectra for each star with Vrot-macro ranging 
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V rot- m »cro [kmS- 1 ] 

1 2340 1 234 




l i l i l i l i l 

6301.4 6301.6 6302.4 6302.6 

Wavelength [A] 



Fig. 4. Determination of v ro t-macro for Hip 75181. The two Fe i 
lines at 6301.501 A (left panels) and 6302.494 A (right panels) 
have been used. The top panels show the un-normalized x 2 ~ 
values for different values of v ro t-macn> m steps of 0.5kms _1 
and a spline function connecting them. The v ro t- ma cro for 
which the ^-functions are minimized are also indicated. 
The middle panels show the observed spectra (solid circles) 
and synthetic spectra for five different values of v ro t_ m acro 
(0.5, 1, 1.5,2, 2.5 kms~'). The bottom panels show the differ- 
ence between the observed and the synthetic spectra in the mid- 
dle panels. Up-arrows indicate where the residuals converge, 
and down-arrows where they intersect. 

from 0.5 to 5kms~' in steps of 0.5kms _1 . For each v lot -macro 
we estimated the (un-normalized)^ 2 -value (i.e. the added sum 
of squares of the difference between observed and synthetic 
spectra in a certain wavelength interval). 

When satisfactory log ^/-values have been reached the in- 
tersection points (indicated by down-arrows in Fig. [4ji of the 
residuals between observed and synthetic spectra should lie on 
a vertical position equivalent to the level farther from the line 
core (i.e. where one approaches the continuum) and the residu- 
als for the different v ro t-macio converge (indicated by up-arrows 
in Fig. 0). 

Next, the ^-functions between observed and synthetic 
spectra were minimized in order to find the v rot - maC ro ( see t°P 
panels in Fig.@}, using the new log gf- values. We adopted the 
average v rot _ macro as derived from the two Fe lines whenever 
possible (sometimes one of the two was affected by telluric 
lines). These values are listed in column 7 in Tabled 

4.3. Solar analysis 

Our solar analysis is based on integrated solar light spectra 
that were obtained by observing the Moon, Jupiter's moon 
Ganymede, and the afternoon sky with CES. The solar UVES 
spectrum is the Moon spectrum that can be obtained from 



Table 3. Our derived Solar oxygen abundances, given on the 
absolute abundance scale (e(O) = logNo/Nn + 12.0). All abun- 
dances have been determined under the assumption of LTE. 
NLTE corrections are done according to the prescription in 
Gratton et al. ( 2000i. In the third column we give the O i triplet 
abundances corrected for NLTE. The two last columns give the 
differences between abundances from this work and the stan- 
dard solar photospheric value, e(O)© = 8.83, from Grevesse & 
Sauval JT9^8l . 



Wavelength 


e(O)© 


e(O)© 


Difference 


(A) 


LTE 


NLTE-corr 


LTE 


NLTE 


6300.304 


8.71 




-0.12 




6363.776 


9.06 




+0.23 




7771.944 


8.83 


8.78 


±0.00 


-0.05 


7774.166 


8.88 


8.83 


+0.05 


±0.00 


7775.388 


8.82 


8.77 


-0.01 


-0.06 



ESO's web-pages 5 and the solar FEROS spectra were obtained 
by observing the afternoon sky. The atmospheric parameters 
for our solar model are listed in Table [l] (see also Bensby et 
al l2003l for further discussion). 

Table lists the abundances we derive for the Sun from 
the different lines using the methods that are described in more 
detail in Sects. |4!4||43] and l4~6l below. For the Oi triplet we 
also give NLTE corrected abundances (see Sect. I4.6I >. There 
are quite large differences between the abundances that the 
different indicators give. The [Oi] lines give e(0)63oo = 8.71 
and e(0)6363 =9.06 while the Oi triplet lines give abundances 
slightly higher than 8.8 dex if no NLTE correction is applied 
and slightly lower than 8.8 dex if corrected for NLTE effects 
according to the prescription of Gratton et al. {1999 1. That the 
abundances from the Oi triplet do not agree with the abun- 
dances from the [O i] lines is not unexpected due to the NLTE 
effects that the triplet lines are known to suffer from (e.g. 
Kiselman fT99H . 

Lambert (1978 1 noted that there is a possible weak CN 
blend in the [Oi]6363 line and estimated its equivalent width 
to be ~ 0.1 8 mA in the Sun. Reducing our measured equiva- 
lent width by this amount results in an oxygen abundance that 
is 0.03 dex lower, i.e. there is still almost a 0.3 dex difference 
between abundances from the two forbidden lines. In order to 
bring down the abundance from the [O i]6363 to what we find for 
the [O i]6300 line we would have to reduce the equivalent width 
from 2.3 mA to 1 .0 mA, i.e. the CN blend would have a strength 
of 1.3 mA in the Sun (if this is the only blend). Another possi- 
ble reason for the inconsistency between the abundances from 
the two forbidden lines is erroneous logg/-values. However, 
the log values for these lines come from the same calcula- 
tions (see Sect. [3}- If they are erroneous it is therefore most 
likely that both lines show deviations in the same direction and 
of equal size. Further, the [Oi]6363 line is located in the wing 
of a Ca i auto-ionization line (Mitchell & Mohler 1965 1 which 

5 http : //www . eso . org/ observing/ df o/quali ty/UVES/ 
pipeline/solar_spectrum.html 
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0.4 



0.2 




-0.6 -0.4 -0.2 0.2 0.4 
[Fe/H] 

Fig. 5. The observed Ni trend for our stellar sample (see 
Bensby et al. 2003 1. The six additional stars have been included 
as well (Bensby et al., in prep.). Thin and thick disk stars are 
marked by open and filled circles, respectively. 



might affect the strength, and thereby the derived oxygen abun- 
dance, from the [Oi]6363 line (see Sect. 14. 5> . 

The standard solar photospheric oxygen abundance is 
e(O)© = 8.83 + 0.06 (Grevesse & Sauval [13381 However, de- 
pending on which spectral line that is used likely values range 
from ~ 8.6 dex to ~ 8.9 dex using one-dimensional (1-D) mod- 
els (Asplund 2003 1. The new 3-D models give better agreement 
between the derived Solar abundance from different indicators 
(Asplund 2003 1. It should be noted that our solar oxygen abun- 
dance derived from the [Oi]63oo line, using standard 1-D LTE 



models, agrees nicely with what Allende Prieto et al. {2001) 
found using 3-D models (e(O)0 = 8.69). 

Since we want to compare the oxygen abundance trends to 
those oxygen trends that can be found in other works we nor- 
malize our abundances individually, line by line, to the solar 
values that we derived. This is done by correcting the abun- 
dances for all stars by subtracting the amounts that are given in 
the two final columns in Table [3] This means that the Sun will 
have [O/H] = independent of which line that is considered. It 
is these corrected oxygen abundances that we give in Tabled 

4.4. Abundances from the [O i] line at 6300 A 

The oxygen abundance from the [Oi]6300 line has been deter- 
mined for 63 of the 72 stars. In 60 cases we used CES spectra, 
in 2 cases FEROS spectra, and in 1 case a UVES spectrum, 
as indicated in Table [2 For 7 stars we were unable to use the 
[Oi]63oo line due to strong interference from the telluric emis- 
sion line (see Sect. 12.31 . In two other cases where abundances 
from this line are lacking we did not have CES spectra and the 
lower quality of the FEROS spectra did not allow any abun- 
dance determination from this line. 

The synthesis of the [Oi]63oo line was done taking the 
blending Nil lines into account. Figure |5] shows the observed 
Ni trend for the stellar sample. It is obvious that [Ni/H] is not 
varying in lockstep with [Fe/H] which is commonly assumed 
in the modeling of the [Oi]63oo line (e.g Nissen et al. l2002t . At 
[Fe/H] > solar-scaled Ni abundances will underestimate the 
Ni contribution to the line and consequently overestimate the 
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Fig. 6. Synthesis of the [Oi]63oo line for three stars. In each 
plot the observed spectrum is marked by solid circles and the 
synthetic spectra by lines: the [Oi] contribution by a dashed 
line; the contribution from the Ni lines by a dotted line; and the 
joint [Oi] - Ni contribution by a solid line. The solar spectrum 
is in this case a Moon spectrum. Abundances as indicated. 



final oxygen abundance. However, at [Fe/H] < the [Ni/Fe] as 
assumption might be reasonably correct as the Ni contribution 
to the total line strength decreases rapidly below [Fe/H] = 0. As 
can be seen in Fig.|(J the Ni i contribution to the combined line 
profile is small but not negligible for Hip 14086 which is a star 
that have [Fe/H] = -0.59. For the Sun the Nil contribution to 
the total absorption is slightly less than 50% (Fig.|5p) and for 
Hip 95447 that has [Fe/H] = 0.37 the Ni i contribution exceeds 
50% of the combined [O i]-Ni i line profile (Fig. [(J). 

To determine the oxygen abundance we first created a set 
of synthetic spectra representative of different oxygen abun- 
dances (in steps of 0.03 dex), using the v rot _ macro determined 
in Sect. 14.21 Each synthetic spectrum was compared to the ob- 
served spectrum and a corresponding (un-normalized)^ 2 -value 
was calculated. In this way we got a ^-function versus oxy- 
gen abundance that was minimized to find the best value of the 
oxygen abundance for each star. The corresponding abundance 
for the Nil blends were kept at a fixed value and was gath- 
ered from Bensby et al. (2003 1. Tables^and|4]lists the derived 
abundances. 



4.5. Abundances from the [O i] line at 6363 A 

The [Oi]6363 line is approximately 50% weaker than the 
[Oi]63oo line and is located in the right wing of a wide Cai 
auto-ionization line (Mitchell & Mohler 1965). This makes 
it difficult to synthesize this line well in our FEROS spectra. 
Instead we determined abundances from equivalent width mea- 
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Fig. 7. Difference between the oxygen abundances from the 
[Oi]6300 line and the [Oi]6363 line versus [Fe/H], logg, and T e s. 
Thin and thick disk stars are marked by open and filled circles, 
respectively. 

surements. That the continuum is depressed is then of small 
concern since the local continuum can be set at every measure- 
ment. However, depending on at which depth in the stellar at- 
mosphere this Ca i auto-ionization line forms its effect on the 
equivalent width of the [Oi]6363 line will vary. If it is formed 
below the [Oi]6363 line it will only contribute to the contin- 
uum flux, while if it forms higher up in the atmosphere than the 
[Oi]6363 line it will increase the strength of the [Oi]6363 line. If 
this latter is the case and if the estimation of the strength of the 
CN blend by Lambert ( 1978 ) is correct (and also assuming that 
CN is the only blend), this implies that the Ca i auto-ionization 
line contributes ~ 1.1 mA to the [Oi]6363 line in the Sun. 

Both [C/H] and [N/H] evolves in lockstep with [Fe/H] at 
high metallicities (Andersson & Edvardsson 1994, Gustafsson 
et al. [13391 Clegg et al. [TWA Shi et al. I5UU21 . From this it 
does not necessarily follow that [CN/H] evolves like [Fe/H]. 
The amount of CN present in a stellar atmosphere depends 
on its dissociation energy and the temperature in the stellar 
atmosphere. The dissociation energy for CN is, however, not 
well known (see discussion in de Laverny & Gustafsson 1998 1. 
But, as can be seen in Figs. and c, there is no trends in 
[O63oo/06363] with either [Fe/H] or T e g which possibly indi- 
cates that [CN/H] evolves with [Fe/H] in our stars. 

Since we have neglected the CN blend (and other poten- 
tial effects) in the abundance determination there is a poten- 
tial risk that the oxygen abundances from the [Oi]6363 line are 
too high. How it will affect the [O/Fe] vs [Fe/H] trends de- 
pends on how the CN abundance evolves with [Fe/H]. It seems 
that [CN/Fe] * for all [Fe/H] (if CN is the only blend) since 
the (normalized) oxygen abundances are the same for both of 
the forbidden lines (see Figs.^t-c). The effect on the derived 



[O/Fe] trends will therefore only be a systematic vertical shift 
for all stars since we neglected the contribution from the blend 
in our Solar analysis as well. 

Due to the lower resolution and the lower signal-to-noise 
of the FEROS spectra compared to the CES spectra we could 
only measure the equivalent width of the [Oi]6363 line in 32 
stars. A UVES spectrum was used in one case (see Table 
The measured equivalent widths for the [Oi]6363 line are listed 
in Tablegl 

4.6. Abundances from the O i triplet at 7774 A 

Abundances were determined through equivalent 
width measurements of the permitted Oi lines at 
7771.95/7774.17/7775.39A in our FEROS spectra for 69 
stars. The UVES spectra for the three additional stars could 
not be used since the triplet lines are not present due to the 
setting of the spectrograph. 

The triplet lines are known to give discrepant abundances 
under the assumption of LTE compared to the [Oi]63oo line 
(e.g. Eriksson & Toft \W79l Kiselman ITWTl [T^l Gratton 
et al. 119991 . An LTE analysis will always give abundances 
that are overestimated. In our analysis we correct the derived 
abundances for NLTE effects according to the prescription of 
Gratton et al. (I1999> . These non-LTE corrections are based 
on comparisons between abundances obtained with statistical 
equilibrium calculations and abundances provided by an LTE 
analysis with the Uppsala MARCS models. 

Table 0] lists our measured equivalent widths as well as 
abundances for individual lines and NLTE corrections from 
Gratton et al. ( 1999 1. Table^gives the mean abundances (both 
LTE and NLTE) from the three lines. 

In Figs, [jj-c we compare the abundances, not cor- 
rected for NLTE effects (but normalized to the Sun), from 
the triplet lines to those from the [Oi]63oo line for the 60 
stars that have both the [Oi]63oo an d the triplet lines ana- 
lyzed (compare Table [Q. In summary, abundances from the 
triplet give higher abundances than what the forbidden line 
gives. The average difference and the standard deviation are 
<[0 6 30o/0 7 774]> = -0.049 + 0.084. The largest differences are 
seen at the highest metallicities ([Fe/H] >0) and highest tem- 
peratures (r e ff >6100K) and, possibly, for logg >4.4. For the 
NLTE corrected abundances the average difference and the 
standard deviation is lower ([O6300/O7774]) = -0.015 + 0.077 
(see Figs. [HJi — f). The trends with [Fe/H] and \ogg are still 
present but somewhat less prominent and the trend with T e g 
has changed significantly. The triplet lines now give increas- 
ingly lower abundances when going to higher T e ff. The NLTE 
corrections of Gratton et al. J199 9 > therefore could be some- 
what overestimated for the parameter space that our stars span. 

A weighted linear regression 6 with log T e g, log g, and 
[Fe/H] as the dependent variables and the difference between 
the oxygen abundances from the [Oi]6300 line and the triplet 



6 For a description of the algorithm see the manual entry for 
the"robustfit" function in the Matlab Statistics Toolbox User's Guide 
(www . mathworks . com) 
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Fig. 8. a)-c) Difference between the oxygen abundances from the [Oi]63oo line an d the Oi triplet at 7774 A versus [Fe/H], 
logg, and r e ff. No NLTE corrections have been applied to the abundances from the triplet. d)-f) Difference between the oxygen 
abundances from the [O i]63oo line and the NLTE corrected abundances (according to the prescription in Gratton et al. 1999 1 from 
the O i triplet at 7774 A versus [Fe/H], log g, and T e ff . Thin and thick disk stars are marked by open and filled circles, respectively. 

Table 4. Oxygen abundances for the 72 program stars. Column one gives the Hipparcos number, columns 2-6 the derived 
oxygen abundances for the different indicators relative to the photospheric oxygen abundance (e(O)0 = 8.83) from Grevesse & 
Sauval ( 1998 1. Columns 7-11 give the abundances relative to our solar abundances, i.e. [0/H]q =0 for all indicators. Columns 
12- 14 give the abundances for the triplet that have been corrected for NLTE effects according the the prescription in Gratton 
et al. (120001 . and normalized to the new NLTE corrected abundance for the Sun. Column 15 gives the NLTE correction terms 
from Gratton et al. (2000 1 (GCEG) that were applied to the abundances in columns 12- 14. Column 15 gives the empirical NLTE 
correction we determined (Eqs.^andllJ. Columns 16-19 give the measured equivalent widths for the [Oi]6363 and the triplet 
lines. The full table is available in electronic form at the CDS via anonymous ftp to cdsarc .u-strasbg . fr (130.79.128.5) or 
via cdsweb.u-strasbg.fr/cats/ J .A+A.all.htx 
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lines ([O6300/O7774]) as the independent gives: 



5. Errors in resulting abundances 



1 7774 



= [(WO7774] = " 0.46 (±0.02) 



- 0.89 (±0.37) -log(T eS /T eSB ) 
+ 0.10 (±0.04) -logte/so) 

- 0.10 (±0.03)- [Fe/H]. (1) 

This empirical correction can then be applied to the abundances 
from the 1 triplet: 



[0/H]n™e - [O/HJltE + A NLTE- 



(2) 



The derived relationship is valid for stars with: 
5000 < r eff < 6500, 3.5 < logg < 4.5, and -1 < [Fe/H] < 0.5. 



The random (internal) errors in the derived oxygen abundances 
can have several causes. Effective temperature, surface grav- 
ity, microturbulence, measured equivalent widths, and atomic 
data are the main contributors. In estimating the total error we 
assume that the errors arising from each of the sources are inde- 
pendent. The error due to each parameter is estimated by vary- 
ing the parameter in question. 

As in Bensby et al. (2003 1 we estimate the effects of the 
uncertainties in these parameters on the derived oxygen abun- 
dances by increasing T e s by 70 K, log g by 0. 1 dex, micro- 
turbulence (£ t ) by 0.15kms _1 , [Fe/H] by 0.1 dex, correction 
factor (dye) to the classic Unsold by 50%, and the measured 
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Table 5. Estimates of the effects on the derived abundances 
due to internal (random) errors for three stars. When cal- 
culating AW A / y/N for the Fe abundances we have assumed 
AW A = 5% for Hip 88622 and Hip 82588 and AW A = 10 % for 
Hip 103682 (see Bensby et al. 2003 1. For the [Oi] 6363 abun- 
dances we adopted AW A = 10% and for the Oi triplet abun- 
dances AW A - 5% for all three stars. The total random uncer- 
tainties (cr ran( j) were calculated assuming the individual errors 
to be uncorrected. The final line gives the average of the total 
random uncertainties for the three stars. 
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equivalent widths by 5 - 10 %. This uncertainty in the measured 
equivalent widths correspond to an uncertainty in the derived 
abundances of 0.02- 0.04 dex. 

If the abundance is derived from spectral lines the un- 
certainty should be reduced by a factor N 1/2 to give the formal 
error in the mean. For the [Oi]6363 line we adopt the higher 
value (0.04 dex) since the line is very weak and difficult to mea- 
sure. For the Oi triplet abundance we adopt the lower value 
(0.02/ V3 dex) since the three lines are strong and easy to mea- 
sure. 

The uncertainty in the abundance from the [Oi]63oo line is 
also dependent on the accuracy of the derived Ni abundance. 
From the discussion in Sect. l4.4l it is also clear that stars at high 
[Fe/H] are more affected by this than those at low [Fe/H]. To 
estimate these uncertainties we vary the Ni abundance accord- 
ing to the observed formal error in the mean of the Ni abun- 
dance (see Bensby et al. 2003) and then derive new oxygen 
abundances from line synthesis. 

These error estimates are exemplified for three stars and are 
listed in Table |5] It should be noted that the error estimates for 



the damping constants only apply to those lines marked by a 
"U" in Table|2] The total internal (random) errors are typically 
0.06 dex in [0/H] 630 o and [0/H] 636 3, and 0.08 dex in [0/H] 77 74- 
We note that the oxygen trends in Fig. and b are very 
tight and that the two populations separate very well. That this 
separation should be so distinct and the internal errors as large 
as the ones estimated in this section seems unlikely. We there- 
fore conclude that our error estimates in Table |5] are upper lim- 
its to the internal errors. 



6. [O/Fe] vs [Fe/H] trends 

Our oxygen trends derived from the [Oi]63oo line are shown 
in Fig. The thin and thick disk trends are well separated 
and distinct. At [Fe/H] < -0.4 the [O/Fe] trend is nearly flat 
for the thick disk stars, showing an oxygen overabundance of 
[O/Fe] as 0.4. Then, at [Fe/H] a -0.4, there is a downward trend 
towards solar values. This trend is in good agreement with what 
we found for other a-elements (Mg, Si, Ca, and Ti) (see Bensby 
et al. 2003 Feltzing et al. l2003> . As we already noted then, this 
turn-over or "knee" is most likely an indication of contribu- 
tion from SN la to the chemical enrichment of the interstellar 
medium. 

The thin disk [O/Fe] trend is very well defined with a 
small internal scatter. It shows an overabundance of oxygen of 
[O/Fe] a 0.3 at [Fe/H] su -0.7 and continues linearly down to an 
underabundance of [O/Fe] w -0.2 at [Fe/H] w 0.4. Such a trend 
has previously only been seen by Castro et al. (|T997 1 who ana- 
lyzed the [O i]6300 line in seven stars with +0. 1 < [Fe/H] < +0.5, 
and by Chen et al. (2003 1 who analyzed the Oi triplet in 15 
stars with -0.1 < [Fe/H] < +0.5. However, only seven stars in 
the Chen et al. (2003) study had [Fe/H] > +0.1. 

Compared to the abundance trends obtained from the 
[Oi]6 3 oo line, Fig.|9^, the trends based on the [Oi]6 3 63 line is 
practically identical even though the number of stars is consid- 
erably smaller. Since the [O i]6363 line of a star is not affected by 
its Ni abundance this provides an independent test of our analy- 
sis of the [Oi]63oo line. We therefore believe that the downward 
trend in [O/Fe] that we see at [Fe/H] > as well as the separa- 
tion of the thin and thick disk trends at [Fe/H] < are real. 

In Fig. [9]: we show the abundance trends derived from the 
O i triplet that have not been corrected for NLTE effects. At 
[Fe/H] < the trends for the thin and the thick disks are essen- 
tially the same as we saw from the forbidden lines (Figs. [9^ 
and b) suggesting that the NLTE effects might not be so severe 
at these metallicities (see also Fig. |8^). However, the scatter 
is significantly larger, especially for the thin disk. At super- 
solar [Fe/H] there is a large scatter and no clear sign of the 
downward trend that we see in the [O i] lines. The [O/Fe] trend 
based on the NLTE corrected (according to the prescription by 
Gratton et al. 19999 Oi abundances is shown in Fig.|9jl. While 
not much have changed at [Fe/H] < 0, except for a somewhat 
more well-defined thin disk trend, the scatter is considerably 
reduced at [Fe/H] > and with a marginal downward trend in 
[O/Fe] . However, the forbidden lines always give superiorly de- 
fined trends. 
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Fig. 9. The oxygen trends for our stellar sample as it appears when using; a) the [Oi]63oo line, b) the [Oi]6363 line, c) the Oi 



triplet at 7774 A, and d) the O i triplet at 7774 A corrected for NLTE effects according to the prescription in Gratton et al. ( 1999 1 
(but see discussion in Sect. 14. 6> . Thin and thick disk stars are marked by open and filled circles, respectively. The error-bars in 
the upper right corner are the total random errors (see Sect.|5]and Table[5}. 



7. Discussion 

In Fig.^|we show the oxygen trends, both the common [O/Fe] 
vs. [Fe/H] trend and the nowadays often used [Fe/H] vs. [O/H] 
trend, with all 72 stars in the sample included. The abundances 
from the [Oi]63oo line have been used whenever available. If 
not, our first choice are the abundances from the [Oi]6363 line 
and second the NLTE corrected (according to the prescription 
of Gratton et al. 1999 1 abundances from the O i triplet. 

7.1. Interpretation of the oxygen trends in the thin and 
thick disks 

7.1 .1 . Strengthening of the trends 

Nissen et al. (2002i analyzed oxygen in 41 stars. Out of those, 
3 1 stars have kinematic data available. For these we have cal- 
culated their TD/D and TD/H ratios and selected those with 
TD/D > 1 and TD/H > 1 as thick disk stars and those with 
TD/D < 1 and TD/H > 1 as thin disk stars. Four stars that have 
TD/D > 1 and TD/H < 1 were consequently interpreted as halo 
stars. Apart from the 4 halo stars this left us with 10 thick disk 
and 17 thin disk stars from Nissen et al. (12002 >. In Fig.llllwe 
add these stars to our own data set. 

The thick and thin disk trends seen in Fig. HOh are fur- 
ther strengthened by the inclusion of the stars from Nissen 
et al. (2002 1 (see Fig. II 1>. Especially so the flat trend in the 
thick disk up till [Fe/H] = -0.5 and the break in [O/Fe] at 
[Fe/H] « -0.4. The thin disk trend at [Fe/H] < is nicely sup- 
ported too. Intriguingly, the most metal-rich halo star is low 
in [O/Fe] compared to thick disk stars at the same metallic - 



ity, which was also found by Nissen & Schuster d!997> . They 
speculated that these anomalous halo stars can have been ac- 
creted from dwarf galaxies that have a different chemical his- 
tories than the Galactic disks. 

For the metal-rich disk ([Fe/H] > 0) the Nissen et al. d2"0"0"2l 
stars have higher oxygen abundances than ours and level out 
to a constant value of [O/Fe]~0. The Nissen et al. (2002i 
stars include a re-analysis of the Nissen & Edvardsson ( 1992 1 
data that did not account for the Ni blend in the [Oi]63oo line. 
Nissen et al. ( 2002 1 calculated the equivalent width of the Ni 
blend assuming solar-scaled Ni abundances (i.e. [Ni/Fe] = 0) 
for all metallicities (which is not consistent with what we see 
in our Ni trends, Fig.|5] but is consistent to earlier studies e.g. 
Edvardsson et al. 1993 1. Then they subtracted this calculated 
Ni contribution from the observed equivalent width of the [O i] 
line and then re-calculated the oxygen abundances. Apart from 
an absolute shift in the oxygen abundances, their new results 
showed essentially no differences compared to the abundance 
trends in Nissen & Edvardsson ( 1992 1. This is not surprising 
since the assumption that [Ni/Fe] = for all [Fe/H] mainly re- 
sults in a systematic shift of all abundances by a certain amount 
and will not affect the behaviour of the [O/Fe] trend. Hence the 
Nissen et al. ( 120021 [O/Fe] trend still level out at [Fe/H] > 0. 

7.1 .2. Evolution of the thick disk 

In Bensby et al. (2003 1 and Feltzing et al. (2003 1 we discuss 
different formation scenarios for the thick disk. Based on our 
elemental abundance trends in the thin and thick disks for the 
other (T-elements (Mg, Si, Ca, and Ti) we concluded that it is 
very likely that the thin and thick disks formed at different 
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Fig. 10. The oxygen trend, as it appears when all 72 stars are included, a) shows the [O/Fe] vs. [Fe/H] trend and b) the [Fe/O] 
vs. [O/H] trend. Abundances derived from the [Oi]6300 line are marked by circles, abundances from the [Oi]6363 line by squares, 
and NLTE corrected (according to the prescription in Gratton et al. 2000 1 abundances from the O i triplet at 7774 A by triangles. 
Thin and thick disk stars are marked by open and filled symbols, respectively. 



epochs. This conclusion is further supported by the different 
mean stellar ages that we found from isochrones for the thin 
and thick disk stars. The average ages are (including all stars 
in this study) 4.9 + 2.8 and 10.8 + 4.3 for our thin and thick 
disk stellar samples, respectively. However, our results could 
be accommodated both in a dissipational collapse scenario (ei- 
ther fast or slow) or in a merger/interacting scenario for the 
formation of the thick disk. However, evidence from extra- 
galactic studies of edge-on spiral galaxies (e.g. Reshetnikov 
& Combes [T9971 Schwarzkopf & Dettmar 2000] Dalcanton & 
Bernstein 2002 1 appear to indicate that thick disks are more 
common in galaxies that are or have experienced mergers or 
interaction. Thus, at the moment the most plausible forma- 
tion scenario for the thick disk is a merger event (e.g. Quinn 
et al. I1993i or an interaction with a companion galaxy (e.g. 
Kroupa l2002l 

The thick disk shows the signatures of enrichment from 
SN la to the interstellar medium. This is clearly seen in plots 
where [O/Fe] is plotted versus [Fe/H] (Fig. llOh and ll 111 . The 
"knee" in [O/Fe] at [Fe/H] « -0.4 indicates when the SN la rate 



peaks and thereby also the peak in the enrichment of Fe from 
these events. That the "knee" is present at all also indicates that 
the star formation in the early thick disk has been vigorous. 
This can also be seen in Fig.llOb that shows how [Fe/O] runs 
with [O/H] as the uprising trend in [Fe/O] at [O/H] ~ 0. 

7.1 .3. Evolution of the thin disk 

The Galactic thin disk has not had such an intense star forma- 
tion history as the thick disk. The shallow decline in the [O/Fe] 
trend from an oxygen overabundance of [O/Fe] ~ 0.2-0.3 at 
[Fe/H] — 0.8 to an oxygen underabundance of [O/Fe] w -0.2 
at [Fe/H] ~ 0.4 indicates a more continuous star formation with 
no fast initial enrichment from SNII. Instead the observed 
[O/Fe] trend favours a chemical evolution in the thin disk where 
both SN la and SN II in a steady rate contribute to the enrich- 
ment of the interstellar medium and thereby producing a de- 
cline in [O/Fe] with [Fe/H] (see Fig.fTUk). 

How can we explain the trend observed in the thin disk? 
As noted in the preceding section the thin disk stars, on aver- 
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Fig. 11. Oxygen trends based on the [Oi]63oo line for our stars 
and from the study by Nissen et al. (2002i. The Nissen et 
al. (2002 1 thin and thick disk stars are marked by open and 
filled triangles, respectively, and their halo stars by asterisks. 
Our thin and thick disk stars are marked by open and filled cir- 
cles, respectively. 



age, are younger than the thick disk stars. However, the thin 
disk stars appear to span the same metallicity range (below so- 
lar metallicity) as the thick disk stars do. A possible scenario 
for the star formation in the thin disk would be that once star 
formation in the thick disk was stopped (truncated?) there is a 
pause in the star formation. During this time there is in-falling 
fresh gas that accumulates in the Galactic plane, forming a new 
thin disk. Also the remaining gas from the thick disk will settle 
down onto the new disk. Once enough material is collected star 
formation is restarted in the new thin disk. The gas, though, has 
been diluted by the metal-poor in-falling gas. This means that 
the first stars to form in the thin disk will have lower metallic - 
ities than the last stars that formed in the thick disk. Although 
this scenario nicely explains the abundance trends in the thin 
disk it remains to work out the details and see if they fit other 
observational tests such as G dwarf and metallicity distribu- 
tions. 

Prior to this study the models of Galactic chemical evo- 
lution did not match the observed [O/Fe] trend at [Fe/H] > 
since the models predict a downward trend (see e.g. Chiappini 
et al. 2003 for the most recent models) while the observed 
trend leveled out (see e.g. Nissen et al. 2002 and discussion 
in Sect. 17. 1.11 . Even though the models of Galactic chemical 
evolution rarely are evolved beyond [Fe/H] as +0.2 there is no 
indication of a leveling out. 

7.2. Synthesis of a-elements and the choice of 

reference elements in chemical evolution studies 

The models of Galactic chemical evolution predict that 
[Mg/Fe], at [Fe/H] >0, should show a downward trend simi- 
lar to that seen for [O/Fe]. Our observations do not show this 
(Bensby et al. |2003 Instead for a downward trend [Mg/Fe] 
levels out at [Fe/H] ss 0. In Fig. [21 we show [O/Mg] as a func- 
tion of [Mg/H] which further emphasizes the different trends 
for oxygen and Mg. 



Prior to the new models by Meynet & Maeder ( 2003} the 
models by Maeder ( 1992 > predicted a strong dependency on 
the metallicity for the oxygen yields. This could then have 
explained the trend in Fig. El But, since the oxygen yields 
for massive stars now have been shown to be less sensitive to 
the metallicity when including rotation in the stellar models 
(Meynet & Maeder 2003 1 this explanation is no longer viable. 
The stellar evolution models are also well-known to underesti- 
mate the Mg yields for massive stars (e.g. Chiappini et al. 2003 
Chiappini et al. 1199 9 Thomas et al. [1998 » and thereby giving 
too low [Mg/Fe] values at high [Fe/H] in the models of Galactic 
chemical evolution. Currently, this could then be the explana- 
tion for the discrepancy between models of Galactic chemical 
evolution and observations for Mg. However, this conclusion 
must remain tentative and more work on stellar yields is neces- 
sary for both Mg as well as oxygen in order to achieve a clear 
picture. In this our observational data can be used to further 
constrain the possible models. 

In studies of chemical evolution of single or multiple stel- 
lar populations it is desirable to have a reference element that 
has only one, well-understood source. Stellar spectra are rich 
in easily measurable Fe lines. This means that for almost all 
stars Fe abundances are readily available and hence Fe is by far 
the most commonly used reference element in studies of galac- 
tic chemical evolution. However, Fe is produced in both SN II 
and SNIa and therefore the trends are less easy to interpret 
than if we had a reference element that has only one source. 
Such reference elements could be provided by Mg or oxygen. 
However, in light of our own observations, the shortcomings 
of stellar yield models, and the fact that at least Mg might be 
produced, to some extent, in SN la we would be very cautious 
about promoting one specific reference element. It is probably 
best, for now, to consider more than one reference element. 

8. Summary 

We have presented oxygen abundances for 72 nearby 
F and G dwarf stars in the solar neighbourhood with 
-0.9 < [Fe/H] < +0.4. The stellar samples were chosen in or- 
der to investigate two important issues related to the formation 
and chemical evolution of the stellar disks in our Galaxy: 

★ the oxygen trends in the thin and thick disk, 

★ the oxygen trend in the thin disk at super-solar metallicities. 

The stars have been selected so that they either with a high 
likelihood belong to the thick or the thin disk. This selection 
was purely based on kinematics. The majority of the stars in the 
study presented here have previously been studied by Bensby 
et al. (2003 1 and Feltzing et al. (2003 1. In those two papers we 
derived stellar parameters, detailed statistics of the kinematic 
selection, and abundances for a- as well as iron peak-elements. 

Our abundance analysis is based primarily on the forbidden 
oxygen line at 6300 A. This line is blended with a doublet of 
Ni lines arising from two different Ni isotopes. As we have 
spectra of very high S /N and R and have modeled the line and 
its oxygen and Ni components in detail we have been able to 
disclose very well-defined trends for oxygen relative to iron 
and magnesium. We have also analyzed the forbidden line at 
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6363 A and the permitted triplet lines around 7774 A and found 
consistent trends from these lines. 

Our main conclusions from this study are: 

(i) At super-solar [Fe/H] the [O/Fe] trend for the thin disk 
continues linearly down-wards. This is different from the 
other ff-elements which show a leveling out of [a7Fe] at 
[Fe/H]=0. 

(ii) The thick disk is more overabundant in oxygen than the thin 
disk at sub-solar metallicities. The thick disk also shows the 
signatures of chemical enrichment from SN la. 

(iii) Oxygen and magnesium do not evolve in lockstep at super- 
solar metallicities. 

(iv) By comparing oxygen abundances from the permitted in- 
frared O i triplet to those from the forbidden line at 6300 A 
we provide an empirical NLTE-correction relation for the 
abundances from the triplet lines that. This could be used 
e.g. for F and G dwarf star spectra with a S/N that is such 
that only the triplet lines that can be analyzed well, e.g. due 
to the distances of the stars. 

(v) Comparing our abundances from the 6300 A line and 
the triplet lines with and without NLTE corrections from 
Gratton et al. ( 1999 1 we find that their NLTE corrections, 
for the parameter space spanned by our stars, are somewhat 
overestimated. 

Conclusion (i) is a new result and indicates that chemi- 
cal evolution models of the Milky Way are now in concor- 
dance with the observed [O/Fe] trends at [Fe/H] > 0. The dif- 
ferent trend for oxygen compared to the other a-elements at 
[Fe/H]>0 (see Bensby et al. 2003 1 indicates that oxygen is 
only produced in SN II. The leveling out of the [a/Fe] trend 
for other <?-elements could be expected on grounds that they 
also have small contributions from SN la. However, given the 
still large uncertainties in SN yield calculations this conclusion 
remains tentative. 

Conclusion (ii) means that the Galactic thin and the thick 
disks indeed have different chemical histories. This is a confir- 
mation and strengthening of previous findings by us and oth- 
ers (e.g. Prochaska et al. l2000> . We also trace the signature of 



SN la (the "knee") in the [O/Fe] trend in the thick disk that 
we have noticed previously for Ca, Mg, Si, and Ti (Bensby et 
al. 120031 Feltzing et 31. 120031 . 
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